function [f,g]=fun(x)
global eta
f=-eta'*log(x);
g=-eta./x;